Lab9

In [1]:

import gtlparser

import gtlparser

Question 1: Creating an Interactive Map¶

In [2]:

map1 = gtlparser.Map()
map1.add_basemap_gui(position="topright")
map1.add_layer_control()
map1.add_search_control(position="topleft")
map1

map1 = gtlparser.Map() map1.add_basemap_gui(position="topright") map1.add_layer_control() map1.add_search_control(position="topleft") map1

Out[2]:

Question 2: Adding Map Legends¶

In [3]:

import leafmap

map2 = leafmap.Map()

wms_url = "https://services.terrascope.be/wms/v2?"
layer_name = "WORLDCOVER_2021_MAP"

map2.add_wms_layer(
    url=wms_url,
    layers=layer_name,
    name="ESA World Cover 2021",
    format="image/png",
    transparent=True,
    shown=True,
    version="1.3.0",
)

map2.add_legend(builtin_legend="ESA_WorldCover")
map2

import leafmap map2 = leafmap.Map() wms_url = "https://services.terrascope.be/wms/v2?" layer_name = "WORLDCOVER_2021_MAP" map2.add_wms_layer( url=wms_url, layers=layer_name, name="ESA World Cover 2021", format="image/png", transparent=True, shown=True, version="1.3.0", ) map2.add_legend(builtin_legend="ESA_WorldCover") map2

Out[3]:

Question 3: Creating Marker Clusters¶

In [4]:

import leafmap
import geopandas as gpd
import pandas as pd  # geopandas relies on pandas

geojson_url = "https://github.com/opengeos/datasets/releases/download/places/wa_building_centroids.geojson"

gdf = gpd.read_file(geojson_url)

gdf
gdf["longitude"] = gdf.geometry.x
gdf["latitude"] = gdf.geometry.y


map3_1 = leafmap.Map()
map3_1.add_marker_cluster(gdf, x="longitude", y="latitude", layer_name="World cities")
map3_1

import leafmap import geopandas as gpd import pandas as pd # geopandas relies on pandas geojson_url = "https://github.com/opengeos/datasets/releases/download/places/wa_building_centroids.geojson" gdf = gpd.read_file(geojson_url) gdf gdf["longitude"] = gdf.geometry.x gdf["latitude"] = gdf.geometry.y map3_1 = leafmap.Map() map3_1.add_marker_cluster(gdf, x="longitude", y="latitude", layer_name="World cities") map3_1

Out[4]:

In [5]:

map3_2 = leafmap.Map()
map3_2.add_circle_markers_from_xy(
    gdf,
    x="longitude",
    y="latitude",
    radius=5,
    color="red",
    fill_color="yellow",
    fill_opacity=0.8,
    marker_cluster=True,
)
map3_2

map3_2 = leafmap.Map() map3_2.add_circle_markers_from_xy( gdf, x="longitude", y="latitude", radius=5, color="red", fill_color="yellow", fill_opacity=0.8, marker_cluster=True, ) map3_2

Out[5]:

Question 4: Visualizing Vector Data¶

In [6]:

geojson_url = "https://github.com/opengeos/datasets/releases/download/places/wa_overture_buildings.geojson"

map4_1 = leafmap.Map()

geojson_style = {
    "color": "red",
    "weight": 2,
    "fillOpacity": 0,
}

map4_1.add_geojson(geojson_url, style=geojson_style, layer_name="Washington Buildings")

map4_1

geojson_url = "https://github.com/opengeos/datasets/releases/download/places/wa_overture_buildings.geojson" map4_1 = leafmap.Map() geojson_style = { "color": "red", "weight": 2, "fillOpacity": 0, } map4_1.add_geojson(geojson_url, style=geojson_style, layer_name="Washington Buildings") map4_1

Out[6]:

In [7]:

geojson_url = "https://github.com/opengeos/datasets/releases/download/places/las_vegas_roads.geojson"

map4_2 = leafmap.Map()

geojson_style = {
    "color": "red",
    "weight": 2,
}

map4_2.add_geojson(geojson_url, style=geojson_style, layer_name="Las Vegas Roads")

map4_2

geojson_url = "https://github.com/opengeos/datasets/releases/download/places/las_vegas_roads.geojson" map4_2 = leafmap.Map() geojson_style = { "color": "red", "weight": 2, } map4_2.add_geojson(geojson_url, style=geojson_style, layer_name="Las Vegas Roads") map4_2

Out[7]:

In [8]:

map4_3_data = (
    "https://github.com/opengeos/datasets/releases/download/us/us_counties.geojson"
)

map4_3 = leafmap.Map(center=(40, -98), zoom=4)

map4_3.add_data(
    map4_3_data,
    layer_name="US County Area (Census)",
    column="CENSUSAREA",
    scheme="EqualInterval",
    cmap="Blues",
    legend=True,
    legend_title="Census Area (sq mi)",
)

map4_3

map4_3_data = ( "https://github.com/opengeos/datasets/releases/download/us/us_counties.geojson" ) map4_3 = leafmap.Map(center=(40, -98), zoom=4) map4_3.add_data( map4_3_data, layer_name="US County Area (Census)", column="CENSUSAREA", scheme="EqualInterval", cmap="Blues", legend=True, legend_title="Census Area (sq mi)", ) map4_3

Out[8]:

Question 5: Creating a Split Map¶

In [9]:

pre_event_url = (
    "https://github.com/opengeos/datasets/releases/download/raster/Libya-2023-07-01.tif"
)
post_event_url = (
    "https://github.com/opengeos/datasets/releases/download/raster/Libya-2023-09-13.tif"
)

map5 = leafmap.Map()

leafmap.split_map(left_layer=pre_event_url, right_layer=post_event_url)

map5

pre_event_url = ( "https://github.com/opengeos/datasets/releases/download/raster/Libya-2023-07-01.tif" ) post_event_url = ( "https://github.com/opengeos/datasets/releases/download/raster/Libya-2023-09-13.tif" ) map5 = leafmap.Map() leafmap.split_map(left_layer=pre_event_url, right_layer=post_event_url) map5

Out[9]: